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APPENDIX I 
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% Creating a network topology object 

network_topo = topo('init'); % graphically place nodes on screen 

addlink(networkjopo); % graphically connect up nodes 

labelnames(networkjopo); % graphically label nodes 

save network_topo; % save networkjopo for future use 



% Top level procedure to compute paths that optimize use of network capacity 
% inputs: 

% D = traffic demand matrix 
15 % (retrieved from predictions stored in TMS Statistics Repository) 

% network__topo = topology object defining the network topology 
% P = network policy information 

% (matrix of reserved capacity, which indicates links whose use 

% is administratively prohibited or which should not be 

20 % completely allocated) 

% outputs: 

% allocated_paths() = list of paths to set up, to TMS signalling system 



25 C = capacity(network_topo); % retrieve network topology information 

C = C-P; 

saved_C = []; 
saved_SLA= []; 
30 assigned_paths = []; 

round = 0; 

[SLA, S] = create_ordered_sla(D); 
35 F = SLA(1) 

for F = SLA', 

round = round +1; 
saved__C {round} = C; 
40 saved_SLA{round} = F; 

F % display the flow 

W = calc_weights('calcweight2 , ,F,C); 
45 [dist,P] = floyd(W); 
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path = fmdpath(P,F.i,F.j); 

assigned_paths {round}. path = path; 
5 assigned_paths {round}. flow = F; 

if (isempty(path)) 

fprintf(l ? 'no path for flow:\n ( ); F 

else 

10 C = compute_residual_capacity('c - F.bw' 5 path,F,C); 

end 

end 

15 function [W] = calc_weights(func,F,C) 

% function [W] = calc_weights(func,F,C) 

% 

% Compute the weights by calling func on each elt of C 

% func must be of the form double func(Flow F, Capacity_elt c, node i, node j) 
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func = fcnchk(func); 



for i = l:size(C,l) 

for j = l:size(C,2) 
25 W(i j) = feval(func,F ? C(i,j),ij); 

end 

end 

function [w] = calcweight2(F,c,ij) 
30 % function [w] = calcweight2(F,c,ij) 

% basic weight calc 

if(0 = c) 

w = inf; 
35 return; 
end 



% rule out paths that can't hack it 
40 if(F.bw>c) 

w = inf; 
return; 

end 
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w = 1 / (c - F.bw) ; % fill links with most capacity first 

function [C] = compute_residual_capacity(func, path, F, C) 
% function [C] = compute_residual_capacity(func, path, F, C) 
% 

% Update capacity characteristics in C to reflect flow F being 

% allocated along path using function func 

% func should be of the form 

% C_element func(C_element c, Flow F) 



if(length(path)<= 1) 
return; 

end 

func = fcnchk(func 5 'cVF'); 



index = 1 ; 
src = path(index); 
20 index = index + 1 ; 

for index = index: length(path) 
dst = path(index); 

25 C(src,dst) = feval(func ? C(src,dst) ; F); 

src = dst; 

end 

30 function [SLA, S] = create_ordered_sla(D) 

% function [SLA] = create_ordered_sla(D) 
% takes the demand matrix and returns a list of SLAs, 
% SLA of the form [ struct ; struct ; ... ] where struct is [BW, i, j] 
% S of the form [ [BW, i, j] ; [BW, i, j] ; ...] 



s = D; 



fori = l:size(D,l) 
40 forj = l:size(D,2) 

if(D(ij)~=0) 

S = [[D(i,j)ij];S]; 

end 

end 

45 end 
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[Y, I] = sortrows(S,l); 

S = Y(size(Y,l):-l :1,:); % reverse order 

5 

SLA = stmct(W,num2cell(S(:,l));i',num2cell(S(: 5 2)), , j',nurn2cell(S(:,3))); 
return; 

10 function [path] = fmdpath(P,i,j) 

% function [path] = findpath(P) 

% 
% 

15 path = []; 

if Ci — J) 

path - [i]; 

return; 

20 end 

if(0==P(i,j)) 
path=[]; 

else 

25 path = [fmdpath(P,i,P(i j)) j]; 

end 

function [D, P] = floyd(W) 
% function [D, P] = floyd(W) 
30 % given weights Wij, compute min dist Dij between node i to j 

% on shortest path from i to j, j has immeadiate predecessor Pij 

n = size(W,I); 
if(n~=size(W,2)) 
35 error('Input W is not square??! !'); 

end 

D = W; 

40 P = repmat([l:n]',[l n]); 

P = P * ~isinf(W); 
P = P * ~eye(n); 

for k = 1 :n 
45 for i = 1 :n 



37 



for j = l:n 

alt_path = D(i,k) + D(k,j); 
if(D(i,j)>alt_path) 

D(i j) = alt_path; 

P(ij) = P(kj); 

end 



end 

end 

k; 

10 D; 

P; 

end 
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APPENDIX II 

function addlink(TOPO) 
% addlink(TOPO) 
% 

% interactively add links to the TOPO 

update(TOPO); 
c_src = 1; 
c_dst = 2; 
c_bw = 3; 

figure(TOPO.cur_fig) 
while (1) 

fprintf(l,'\n\nHit Button 3 to end...\n\n'); 

% find coords and index i of src 
[xli yli button] = ginput(l); 
if (button = 3) break; end 

d = sqrt((TOP0.1ocs(:,l) - xli). A 2 + (TOP0.1ocs(:,2) - yli). A 2); 
[d,i] = min(d); 

xl = TOPO.locs(U); yl = TOP0.1ocs(i,2); 

% find coords and index j of dst 
[x2iy2i] = ginput(l); 

d = sqrt((TOP0.1ocs(:,l) - x2i). A 2 + (TOP0.1ocs(:,2) - y2i). A 2); 
[do] = min(d); 

x2 = TOP0.1ocs(j,l); y2 = TOP0.1ocs(j,2); 
hold on; 

lh = line([xl x2],[yl y2], 'coloured'); 
cap = input('Enter capacity (in Mbps) > '); 

fprintf(l, 'About to create symetric %d Mbps link from node %d to node %d\n',cap,ij); 

doit = input('Enter Y to confirm, N to reject, and B to change bandwidth (Y)> ','s'); 

if (isempty(doit)) doit = 'Y'; end 

if (doit == 'n' | doit == 'N') 
delete(lh); 
return; 
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end 



if (doit = 'b f | doit = TB') 

buf = sprintf( f Enter capacity from %d to %d (in Mbps) > ',i j); 
cap__i_toJ = input(buf); 

buf = sprintf('Enter capacity from %d to %d (in Mbps) > ' j,i); 
capJ_to_i = input(buf); 



else 

1 0 cap_i_to J = cap; 

capj _to_i = cap; 

end 

%build the link records 
1 5 clear linkab linkba; 

linkab.src = i; 
linkab.dst = j; 
linkab. bw = cap_i_toJ; 
20 linkab.handle = lh; 

linkba.src =j; 
linkba, dst = i; 
linkba.bw = cap J_to_i; 
25 linkba.handle = Ih; 



% now draw the actual link on the map 
delete(lh); 

30 lh - drawlink(TOPO, linkab); 

% now store the link info 

TOPO.links = [TOPO.links ; linkab ; linkba ]; 

TOPO.linkarray = [TOPO.linkarray ; [ i j cap_iJoJ] ; [ j i capjjoj ]]; 



end % of while loop 
assignin( 1 caller',inputname( 1 ) ? TOPO); 



function [C, portmap] = capacity(TOPO) 
45 % [C, portmap] - capacity(TOPO) 
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% portmap maps indices of C to elts of nodes(TOPO) 
% [node dir] where 

% node is index of elt in nodes(TOPO) 

% dir is 1 if data enters here, -1 if data leaves here 

5 

numnodes = length(TOPO.links) * 2; 
C = zeros(numnodes ? numnodes); 

10 curnode = 0; 

portmap = []; 

for i = l:length(TOP0.1inks) 
link = TOPO.links(i); 
curnode = curnode 4- 1 ; 
1 5 portmap(curnode 5 :) = [link.src -1]; 

curnode = curnode + 1 ; 
portmap(curnode ? :) = [link.dst 1]; 

C(curnode-l, curnode) = link.bw; 

20 end 

cjiode = 1 ; 
c_dir = 2; 

25 for i - 1 :length(TOPO.nodes) 

ins = find(portmap(: 5 c_node) = i & portmap(:,c_dir) = 1); 
outs = find(portmap(:,c_node) = i & portmap(:,c_dir) = -1); 

for j = ins 

30 for knouts 

CG,k) = inf; 

end 

end 

end 

35 function [a, b, c] = debug(t) 

update(t); 
fieldnames(t) 

40 

a = tnodes 
b = t.locs 
c = t. links 

function display(TOPO) 
45 % DISPLAY a topo object 
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% a link is a unidirectional, so the value is probably twice what you want 

fprintf('[TOPO object: %d nodes %d links]\n',... 
length(TOPO.nodes),length(TOP0.1inks)); 

function draw(TOPO) 
% draw(topo) 

% 

% draw the topology figure in a new window 

TOPO.cur_fig = figure; 

axis(TOPO.axis); 

axis equal; 

axis manual; 

box on; 

hold on; 

for i = l:length(TOPO.nodes) 

nm = plot(TOPO.nodes{i}.loc(l),TOPO.nodes{i}.loc(2),'ob'); 

TOPO.nodes{i}.mark_handle = nm; 
if (isfield(TOPO.nodes{i},'nameloc')) 

TOPO.nodes{i}.nameloc(3) = text(TOPO.nodes{i}.nameloc(l),.. 

TOPO.nodes{i} .nameloc(2),TOPO.nodes{i} .name) 

end 

end 

% yes, this draws the same link twice, fix it if it matters -dam 11/21 

TOPO.linkarray - []; 

for i = l:length(TOP0.1inks) 

TOP0.1inks(i).handle = drawlink(TOPO,TOP0.1inks(i)); 

TOPO.linkarray = [TOPO.linkarray ; ... 

[ TOP0.1inks(i).src TOP0.1inks(i).dst TOP0.1inks(i).bw]]; 

end 

assignin('caller',inputname( 1 ),TOPO); 
function ex(t) 

t.nodes 

function labelnames(TOPO) 
% function labelnames(TOPO) 
% make it easy to label the nodes 

for i = 1 :length(TOPO.nodes) 

fprintf('Place label for node %d "%s M \ri,i,char(TOPO.nodes{i}.name)); 
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origcolor = get(TOPO.nodes{i}.mark_handle,'color'); 
set(TOPO.nodes{i}.mark_handle;color',[l 0 0]); 

if(isfield(TOPO.nodes{i};nameloc')) 
5 good_x = TOPO.nodes{i},nameloc(l); 

good_jr = TOPO.nodes{i}.nameloc(2); 

end 

while (1) 

10 fprintf( ( Button 1 to (re)place text, Button 3 to accept\n'); 

[x,y,button] = ginput(l); 

if (3 == button) break; end 

if (-isempty(th)) delete(th); end 

th = text(x,y,TOPO.nodes{i}.name); 
1 5 good_x = x; good_y = y; 

end 

TOPO.nodes{i}.nameloc = [good__x ? good_y 5 th]; 
set(TOPO.nodes {i} .mark^handle/color'yOrigcolor); 
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end 



assignin('caller',inputname(l ),TOPO);function names(TOPO) 
% NAMES the list of names of the nodes in the topo 



fprintfCNode^tName^'); 
25 for i = 1 :size(TOPO.names,l) 

Q)rintf( , %d\t\t 0 /os\n , J i,TOPO.names{i}); 

end 

function [node] = nodes(TOPO) 
% function [node] = nodes(TOPO) 
30 % returns a cell array describing nodes in the TOPO 

node = TOPO. nodes; 
function [TOPO] = topo(TOPO) 
%[TOPO]=topo(TOPO) 
35 %% if input TOPO is 'init', create a new topology 

% 

% newtopo = topo( f init f ); 

% 

% else add new nodes to TOPO 
40 % 

% nodes is a array of structs, one per node 

% link is a array of structs, one per link 

% a link is a unidirectional item, so there are probably twice 

% as many links as you'd expect. 

45 
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if (nargin < 1) 

error('topo(TOPO) or topo( M inif) - not enough args f ); 

end 

if (ischar(TOPO) & TOPO = 'init') 
clear TOPO 

TOPO.nodes = []; 
TOPO.links - []; 

TOPO .capacity = []; % now computed as needed 
TOPO. Iocs = []; % internal cache 
TOPO iinkarray = []; % internal cache 



15 f= figure; 

axis([0 75 0 50 ]); 

TOPO.axis = axis; 

TOPO.cur Jig - f; 

axis equal 
20 axis manual 

box on 

hold 



else 

25 end 



figure(TOPO.curJig); 



nodecount = length(TOPO. nodes); 



while (1) 

clear nodeinfo; 

fprintf(l,'\n\nHit Button 3 to stop\n\n'); 
[xybut] = ginput(l); 
35 if (but —3) break; end 

x = floor(x); y = floor(y); 
nm = plot(x,y,'ob r ); 
name = input('Enter name > 7s'); 

40 nodeinfoioc = [ x y]; 

nodeinfo. mark_handle = nm; 
nodeinfo.name = cellstr(name); 
nodecount = nodecount + 1 ; 
TOPO.nodes {nodecount} = nodeinfo; 

45 end 
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if('topo' ~= class(TOPO)) 

TOPO = class(TOPO;topo'); 

end 

5 

if (nargout = 0) 

assignin('caller T ,inputname(l) ; TOPO); 

end 

function lh = drawlink(TOPO, link) 
10 % assumes TOPO.linkarray is already valid, and draws the position of 

% link line based on the number of links already present in linkarray 

c_src = 1 ; 
c__dst ~ 2; 
15 c_bw = 3; 

i = link.src; 
j = link.dst; 

20 xl - TOPO.nodes{i}.loc(l); 

yl =TOPO.nodes{i},loc(2); 

x2 = TOPO.nodes{j}.loc(l); 
y2 = TOPO.nodes{j}loc(2); 

25 

if (isempty(TOPO.linkarray)) 
num Jinks ~- 0; 

else 

numjinks = sum(TOP0.1inkarray(: ? c_src) == i & TOP0.1inkarray(: 5 c_dst) = j); 

30 end 

pattern=[0 1-12-2 3-3]*.3; 

if(abs(xl -x2)>abs(yl -y2)) 
35 delta_x = 0; 

delta_y = pattern(num Jinks +1); 

else 

delta_x = pattern(num Jinks +1); 
delta_y = 0; 

40 end 

lh = line([xl x2] + delta_x, [yl y2] + delta_y, 'color', 'black'); 
function update(TOPO) 

45 
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clear TOPO.locs; 

for i = 1 :length(TOPO.nodes) 

TOP0.1ocs(i,:) = TOPO.nodes{i}.loc 

end 

clear TOPO.linkarray; 

for i = l:length(TOP0.1inks) 

TOPO.linkarray = [TOPO.linkarray ; ... 

[ TOP0.1inks(i).src TOP0.1inks(i).dst TOP0.1inks(i).bw]]; 

end 

% these are here to be cut and pasted into other functions as needed 

% there doesn't seem to be a good way to pass them around in another fashion 

% (using assigning('caller'...) to force their definition sounds like asking 

% for trouble 'cause you'll overwrite another definition of them...) 

c_src = 1 ; 

c_dst = 2; 

c_bw = 3; 

assignin( , caller , ? inputname(l) ; TOPO); 



